Space Charge Influence on the Angle of Conical Spikes 
Developing on a Liquid-Metal Anode 
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The influence of the space charge of ions emitted from the surface of a conical spike on its shape 
has been studied. The problem of the calculation of the spatial distributions of the electric field, ion 
velocity field, and the space charge density near the cone tip has been reduced to the analysis of a 
system of ordinary differential equations. As a result of numerical solution of these equations, the 
criterion of the balance of the capillary and electrostatic forces on the conic surface of a liquid-metal 
anode has been determined. It has allowed us to relate the electrical current flowing through the 
system, the applied potential difference and the cone angle. We have compared the results of our 
calculations with available experimental data concerning emission from the surface of pure liquid 
gallium (Ga), indium (In), tin (Sn), and some liquid alloys, such as Au+Si, Co+Ge, and Au+Ge. 
On the basis of the proposed model, explanations have been given for a number of specific features 
of the emissive behavior of different systems. 

PACS numbers: 41.20.Cv, 52.59.Sa 
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Introduction 

It is known that, as a result of the electrohydrodynamic 
instability development, the surface of a conducting liq- 
uid (liquid metal) takes the shape of a cone in a rather 
strong external electric field [l|-|3|- An enhancement of 
the field near the cone tip provides conditions for the ini- 
tiation of emission processes such as field evaporation of 
ions [1411] ■ Interest in studying the geometry of such 
structures was stimulated in large measure by the devel- 
opment of liquid- metal ion sources (LMIS). Considerable 
progress in the theory of conical spikes started from Tay- 
lor's investigations @, [l3j]. He has shown that, for a cone 
with a half-angle of ax — 49.3°, the surface electrostatic 
pressure Pe depends on the distance from its apex as 
R^ 1 and, hence, can be counterbalanced by the surface 
pressure Pl ~ J? -1 . 

The geometry of ion-emitting cones (Taylor cones) 
was investigated in numerous experimental works [bH - 
[l9l |. These works testify that an increase in the applied 
potential difference is accompanied not only by the ap- 
pearance and increase of the emission current, but also by 
a decrease of the cone half-angle. For small currents, i.e., 
when the space charge influence is negligible, the half- 
angle is close to Taylor's angle olt- The phenomenon of 
a cone sharpening can be interpreted as the system re- 
sponse aimed at the conservation of the balance of the 
pressures Pe and Pr, under the conditions of the screen- 
ing effect of the space charge. The last effect in case of 
an emitter with invariable shape, as is known [2fl l2lj |. 
reduces to the Child-Langmuir law. 

Simple analytic models, which relate the basic param- 



*Electronic address: grey@iep.uran.ru 
^Electronic address: nick@ami.uran.ru 
t Electronic address: olga@ami.uran.ru 



ru 



eters of a problem, play an important role in gaining an 
insight into the physical processes that occur in liquid- 
metal ion sources. Among these models are Mair's the- 
ory [22[ . th e models by Kingham & Swanson [Io|. Mair & 
Forbes tHH], and > certainly, Taylor's model [|. Nev- 
ertheless, the above-listed models do not present a theo- 
retical description of the current dependence of the cone 
angle. In this work, we propose a model, which general- 
izes, on the one hand, Taylor's solution to the case where 
the space charge starts playing an important part, and on 
the other hand, the Child-Langmuir law to the case of an 
emitter of variable (self-adjusting) shape. It is based on 
self-similar solutions for a charged particle flow [IE , 
which were found to be compatible with the Laplace- 
Young stress condition for a liquid conducting cone: see 
our Letter [27]]. Clearly, this model does not attempt to 
describe all the features of operation of liquid-metal ion 
sources. The application of self-similar solutions restricts 
our analysis to precisely conical shape of the emitter. So, 
possible deformations of the emitting cone observed in 
experiments, in particular, the appearance of small jet- 
like protrusion on its vertex remain beyond the scope of 
this paper. The advantage of the proposed model is the 
possibility to describe distributions of the electric field 
potential and of the ion velocity field over the cone and, 
as a consequence, to obtain relations between the cone 
angle, the current, and the applied voltage. 



INITIAL EQUATIONS; SELF-SIMILAR 
REDUCTION 



Let us consider a single-velocity flow of ions evapo- 
rating from the surface of an infinite conical anode in 
the framework of the hydrodynamic description [2514281 . 
Figure [1] shows the geometry of the problem and the no- 
tations used. The mass and charge densities are pro- 
portional to each other for a charged-particle flow, and 
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FIG. 1: The geometry of the problem and the notations: L 
is the characteristic spatial scale, a is the cone half-angle, 
m = 7T — a, is the absolute value of electric field potential 
at the symmetry axis at the distance L from the apex of the 
cone. The solid lines correspond to the equipotential surfaces, 
including the surface of the cone with $ = 0. The dashed lines 
correspond to the streamlines of the ion flow. 



vorticity of the flow is conserved in a potential electric 
field (E = — V$). We assume the initial ion velocity to 
be equal to zero. This approximation can be considered 
to be well founded in view of the extremely high values 
of electric field strength near the cone tip. Then the vor- 
ticity of the flow is zero close to the anode surface, and, 
hence, the ion flow is potential. In this case, the flow 
velocity, along with the electric field, can be described 
with the help of a scalar function, the velocity poten- 
tial: v = V^F Then the original set of equations can be 
written as follows: 



V 2 $ = -qN/e , 

m(V*) 2 /2 = -q$, 
V = 0. 



(1) 
(2) 
(3) 



The first equation is the Poisson equation for the electric 
field potential $; here N is the charged particle density, q 
is the particle charge, and Eq is the vacuum permittivity. 
The motion of particles is described by the second equa- 
tion, which gives the energy conservation law for ions in 
an electric field (\I/ is the velocity potential, m is the mass 
of a particle) . The third equation is the continuity equa- 
tion. The Poisson equation ([T]) has to be solved together 
with the condition that the conductor surface is equipo- 
tential. It is convenient to take the electric potential of 
the the anode be equal to zero, $ = 0, which agrees with 
the condition that the initial velocity of particles is zero, 
V* = (see equation (J2j)). 

The equilibrium configuration of the free liquid-metal 
surface is determined by the balance condition for the 



electrostatic and capillary forces (the Laplace- Young 
equation). For a conical surface it takes the form: 



Pt. 



R 



■ cot a 



so 



(V<f>)i =0 = P, 



(4) 



where a is the cone half-angle, a is the surface tension, 
and R is the distance from the cone apex. In the case 
of a high mass flow rate, characteristic for the so-called 
cone-jet mode of electrosprays [2!| |3(| or for sonic sprays 
[3lj ] , the hydrodynamic term also contributes to the pres- 
sure balance condition. However, for the process of field 
evaporation of ions realized in the experiments (T^4l9| . 
this term appears to be 2-3 orders of magnitude less than 
the electrostatic (Pe) and capillary (Pl) terms [3J], and 
we can neglect it. 

In order to pass to dimensionless variables, we intro- 
duce some characteristic spatial scale L, which will be 
taken as the unit of length (see Fig. [T]). Since the equa- 
tions ((T])-© have no their own characteristic spatial size 
(they are invariant with respect to scaling), the value of 
L cannot be defined from the model. Below we will treat 
L as a size of the top part of an infinite model cone; the 
emission current from this part of the cone will be iden- 
tified with the experimental current from a liquid-metal 
ion source. As an external control parameter, we intro- 
duce the absolute value of electric field potential at the 
symmetry axis at the same distance L from the apex of 
the cone, > 0. Then it is convenient to define the 
dimensionless variables by the following way: 



(5) 



r = R/L, ^ = -*/* £ , 

n = NL 2 q/ (s ^l) , </> = * [m/ {2qL 2 <S> L )\ 1/2 , 

and the initial equations (H}-© takes the form: 

V 2 = n, (V^) 2 = cf>, V (nVT/>) = °- ( 6 ) 

Switching to spherical coordinates and taking into con- 
sideration the axial symmetry of the problem, we obtain: 



cot 6 



t>e = n, 



n r ip r 



nglpg 



lp rT 



2 , . 1p66 
-Ipr + — 



cot ( 



■ 1p8 



0. 
(7) 

It can readily be seen that these equations are invariant 
under the scale transformations 



n ■ s 



7-2 



ip if) • s 7 / 2+1 , r — > r ■ s, 



(8) 



where 7 and s are some constants. Consequently, the 
system ([7]) can be reduced by means of the following self- 
similar ansatz (see also [25l |26|): 

= r 7 ^(0), n = r"<- 2 B(6), ip = r 7/2+1 C((9), (9) 
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where A, B, and C are unknown functions of the polar 
angle 0. 

In the dimensionless form, the pressure balance condi- 
tion ((3]) reads 



(Y7A\ 2 C ° ta 

(V0) = O = ~y T) 



V = <i>i 



so 

2Ld 



(10) 



The dimensionless group V plays the role of an external 
control parameter of the problem. One can see from this 
expression that the electric potential has to depend on 
the distance r as </> ~ r 1 / 2 that uniquely determines the 
value of the self-similarity parameter: 7=1/2 0, Il3j . 
Thus, we should apply the substitution 



(/>(r,0)=r L ^A(0), 



= r 1 ' 2 , 
(r,<9) = r- 3/2 B(9), 



(11) 



- r 5/4 



C(0), 



which provides separation of variables in ([7]), The sub- 
stitution (jlll) enables us to reduce the partial differential 
equations (J7J to the following set of second-order ordi- 
nary differential equations for the angle distribution of 
A{0), B{0) and C{0): 

(3/4) A + Agg + Ag COt0 = B, 

(25/16)C 2 + C 2 g = A, (12) 

(15/16)SC + BgCg + BCgg + BCg cot = 0. 

Note that the last equation can be simplified by intro- 
ducing an auxiliary function D(0) = BCg sin#: 



DgCg = -(15/16)1X7. 



(13) 



Solutions of Eqs. (|12p must satisfy a number of con- 
ditions at the symmetry axis (0 = 0) and on the cone 
surface (0 = m = tt — a): 



A(0) = 1, C e (0)=0, A{0 m )=Q. 



(14) 



Note that the condition that the initial velocities of emit- 
ted particles are equal to zero, i.e., the pair of condi- 
tions C(0 m ) = and Cg(0 m ) — 0, follows from the sec- 
ond equation of the set (|12p . Differentiating the same 
equation with respect to 0, we find that the condition 
Ag (0) = is satisfied if the second derivative of the func- 
tion C{0) is finite at the symmetry axis. So, the solution 
of the problem (fT2"]) and (fT4")) is unique for a given value 
of m . 

The additional boundary condition 



An 



= —V 1 Vcot a, 



(15) 



resulting from the pressure balance condition (jlOl) , is not 
used in seeking a solution of Eqs. (fl"2")) . It allows us to 
determine the value of the parameter V that corresponds 
to the angle m . Now, let us rewrite this condition in the 
following form: 



V = Vcata/E(a), E(a) = -Ag{0 m ), 



(16) 



The function E(a) defines the dimensionless value of the 
electric field strength on the cone surface at a distance 
r = 1 from the apex (E(a) = |V</>| r=1 g = g )■ 

The intensity of the ion flux can be characterized by the 
electrical current / (in dimensional form) flowing through 
the sphere of radius L with a center at the apex of the 
cone: 



I = 2nqL 2 / Nt> R sm0d0. 



(17) 



As a corresponding dimensionless quantity, it is conve- 
nient to take the group 



J=I 



(32g 2 e 0J L 3 ( r 3 ) 1 /4' 



(18) 



which does not contain the potential difference $l- It 
can be expressed in terms of the functions B and C by 
following relation: 

J = V 3/2 F(a), E(a) = ( m BC sin0d0. (19) 



The function F(a) defines the particle flux from the top 
part of the cone, < r < 1. 

The self-similar solutions (fTTj) give the ion flux den- 
sity (j ~ n|VV>|) and the electric field intensity (|V</>|) 
proportional to r~ 5 / 4 and, respectively, to r -1 / 2 . This 
corresponds to the power law between the quantities j 
and |V(/>| on the cone surface, 



J 



|V0| 



5/2 



(20) 



It correctly reflects the basic property of the system, 
namely, the nonlinear growth of current density with 
increasing intensity. Nevertheless, the 5/2 power law 
certainly differs from the exponential dependence deter- 
mined by the kinetics of the field evaporation process 
& 

Our solutions would be exact if the actual relation be- 
tween the flux density and the electric field were the same 
as the model relation (|20[) . Note that the self-similar so- 
lutions corresponding to the model law of emission (|2"0")) 
are the only solutions consistent with our main assump- 
tion that the surface of a liquid metal is conical. Any 
distinction of the law of emission from the model law 
(f2T)]) will lead to deviation of the surface from the ideal 
cone. 

However, it is known from the experiments [l4l - fli^ | 
that the surface takes a near-conical shape, and it is 
possible to associate it with a certain cone half-angle a 
(the method of angle measurement is given in the above- 
mentioned papers). That is, the details of the current 
density distribution over the emitter surface (basically, 
ions evaporate from the protrusion growing at the cone 
apex) have a relatively small influence on the balance 
between the electrostatic and capillary forces at the pe- 
riphery of the cone structures. The reason is that the 
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influence of the space charge has an integral character 
(this phenomenon is most conspicuous in planar geom- 
etry, where the screening effect of the space charge in 
principle does not depend on its distribution). As a con- 
sequence, our approach, not claiming to describe LMIS 
operation in detail, is rather applicable for the analysis 
of the integral (averaged) characteristics of LMIS, includ- 
ing the relations between the cone half-angle a, the to- 
tal emission current J, and the applied voltage V. The 
relations between these main model parameters is deter- 
mined by the expressions (fT6"| and p^|) . In the next sec- 
tion, the set of Eqs. (fT2"j) will be solved numerically. It 
will allow us to find the auxiliary functions E(a) and 
F(a), which appeared in fp~6|) and (fT9|). As a conse- 
quence, the dependence of a on J {V plays the role of 
the parameter) and the current-voltage dependence (a is 
the parameter) will be established. 



II. CONSTRUCTION OF SOLUTIONS 
A. Asymptotic expansions 

In order to solve the ordinary differential equations 
(fT2")) and (fT5|) with the conditions (|T4")) numerically, we 
should use asymptotic expansions for the unknown func- 
tions at 9 — ¥ and 9 —> 6 m . This is caused by the 
singular behavior of the functions (or their derivatives) 
at the boundaries 9 = and 9 = 9 m . As will be shown 
below, B — > oo at 9 — > 9 m and Bqqq — >■ oo at 9 — > 0. 

At 9 — > 0, i.e., at the symmetry axis, the functions 
A(6), B(6), C(9), and D{9) can be expanded into the 
series 

A = 1 + a 2 9 2 + a 4 6» 4 + . . . , 

b = b o e + bob^+v + ..., 

C = 4/5 + c 2 9 2 + c 4 9 2 + ..., 
D = d 9 2 +P+dod 2 9 4 +e + .... 



(21) 



Substitution of these expressions into the initial equa- 
tions (H2J), (T3J) yields 



3 1 

«2 = -777, CLi 



(-2 



2 4 ' 2 10 ' 
5 + Vl3 295 - 107^ 



2 4 



C4 



9 ■ 2 



10 



a VL3 + 1 _„ , 275+ 173^13 

p — — Z.o. Pi — 

H 2 ' 1152 

165 + 31^ 



3 . . . , 



5 * * * J 



(22) 



d Q = 2b c 2 , 



d 2 = - 



384 



The expansions (|2Tj) satisfy two first conditions from (|T4|) ; 
they contains a free parameter (the coefficient bo), which 
is determined by the following condition at the cone sur- 
face: A(9 m ) = 0. Note that the divergence of higher 
derivatives of the function B leads to the divergence of 
higher derivatives of the functions A and C. The next 



terms of the expansions (1211) for A and C are of the order 
of0 2 +0. 

In the limit 9 — > 9 m , i.e., on the cone surface, the 
unknown functions can be expanded into power series in 
the parameter x = 9 m — 9: 



A 



,1/2 



E a'iX 



r~i/2 



B 



E W' 2 , 



C = xY,<x 1 / 2 , D = d (l + x 3 / 2 J2d' l x*/ 2 

i=l \ i=l 



(23) 



The first coefficients of these expansions are listed below: 



a 2 = 3 »1, 

, _ 1 b[ 



b' 



d' = — ^i/a^sm 



2 

= ~3 

1 16 



(24) 



Free parameters of these expansions, namely, the coeffi- 
cients a[ and b^, are determined by the boundary condi- 
tions at the symmetry axis. 

Increasing the charge density over the cone leads to 
screening of the electric field at its surface, i.e., the coef- 
ficient ai will approach zero. In the formal limit a[ = 
(the electric field turns to zero at the anode surface), 
another asymptotics is realized: 



A = x 1 / 3 J2a"x i , 



B = x-^ 3 E V(x\ 



C = x 2 / 3 J2 



D = d ' I 1 



x E d'lx 1 

i=l 



(25) 



where 



a i = 1 b i, 

l 1Q V 1) 

b" = — b'l cot 9 ( 
15 

4 = 



a 2 = - b'l cot 9 C 
5 



3 

20 



fy[ cot( 



dl = --b'[ 3 l 2 ^9 c , 
9 

32' 



(26) 



These expressions contain two free parameters: the coef- 
ficient b'l (or tig) an d the angle 9 m = 9 C at which complete 
screening of the external field occurs. 

It should be noted that the function D converges 
rapidly near the cone surface for both forms of expan- 
sion, (|2"3"|) and (|2"5|) . Therefore, using the function D is 
preferable to using the function B in the procedure of 
numerical integration of the problem (|12p - (fT4"]) . 



B. Numerical calculations 

The set of Eqs. (fl"2"|) . (|13p was solved numerically by 
the prediction-correction method of the third order. The 
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FIG. 2: Solutions of the system (|12|l for the cone half-angles 
a = qt (in this case, B = 0) (dotted lines), a = 35° (dashed 
lines), and a = a c (solid lines). 



calculation starts from the asymptotics (|23|) or (|25j) . It 
has been found that the equations admit solutions for 
the angles in the range 9t < 9 m < 9 C . The minimum 
angle value Qt — 130.71° corresponds to the Taylor cone 
(a = olt — 49.29°) and refers to the special case of no 
space charge, i.e., B = 0. For this case, the 9 dependence 
of A is determined by the Legendre function P 1 / 2 (cos9). 
The upper bound of the angle, 9 Cl equals ~ 158.11°. It 
corresponds to the least possible half-angle of the cone 
a = a c ~ 21.89°. This configuration of the surface re- 
lates to the formal limit that the electric field at the cone 
is completely screened by the space charge (this limit 
cannot be achieved because of the finite emissivity of the 
surface). The results of the calculations, corresponding 
to different values of a, are presented in Fig. [5J In view 
of the obvious relations 



Wcf> = r-V 2 {A/2,A e ,0}. 
V7A = r 1 /4{5C/4,C £) ,0}, 



(27) 



the angle dependence of A{9) and C{9) give an idea of 
the distributions of the electric field and the velocity field. 
As an example, Fig. [1] shows the equipotential surfaces 
and streamlines of the flow corresponding to a = 45°. 



C. The electric field strength at the cone surface 

The calculations show that the function E(a) grows 
monotonically from zero to the value ~ 0.975 as the angle 
a increases from a c to olt (see Fig. [3]). Near the limiting 
case (|25l . the estimate E{a) ~ (a — ac) 1 / 2 is valid for 
a — > a c . To verify this dependence, it is possible to 
draw an analogy between the considered problem and 
the problem of a charged-particle flow in a plane vacuum 
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20 30 40 50 

a (deg) 

FIG. 3: The auxiliary functions E(a) and F(a) (solid lines), 
which represent the electric field strength on the cone surface 
at r = 1 and, respectively, the ion flux from the top part of 
the cone (0 < r < 1). The dashed lines correspond to the 
one-dimensional model (Eqs. (|3ip - (|33[) V 



diode, which can be solved analytically. Actually, in the 
one-dimensional case, where all quantities depend only 
on one coordinate z, the equations ([5]) reduce to the set 
of ordinary differential equations 



ip z 2 = 0, (nip z ) z = 0, 



(28) 



which look like the equations (jT2j) after the substitution 
z — > 9. These equations admit the following exact im- 
plicit solution for the electric field potential 4> : 



6j 2 z - E 3 = hjy/j - E 2 ) sj 4j + E Q 



(29) 



where j = nip z is the constant current density, Eo — 
4>z\z=o is the electric field strength at the planar emit- 
ter. Let us also introduce the notation for the electric 
strength at the opposite electrode, E\ = 4> z \ z= fn {h is 
the interelectrode distance) and for the electric potential 
on it, ifih = 4>\ z= h- From the solution (j2T?)) one can obtain 
the relations between the quantities Eq, E±, <f>h, h, and 
3, 



4iv 
6fh - E Q 3 



h — Ei 2 — E 2 , 



Eq ) Ei 



(30) 



By analogy with 9 C , we introduce the interelectrode dis- 
tance h c = 40/j/(3i?i), corresponding to the limit E = 0, 
and then rewrite the expressions (1301) as 



J 



9h 2 



E\ 



En 



Ei + E( t 



(31) 



Note that, in the limit E /Ei —> 0, the first equation 
of (|3ip represents the Child-Langmuir law [2(J [H| for a 
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FIG. 4: Dependence of the cone half-angle a on the dimen- 
sionless emission current J. The inset demonstrates the same 
dependence for relatively small currents. 

plane vacuum diode. It follows from the second equation 
of (f3"Tj) that Eq ps Eiy/l — h/h c in the same limit, i.e., 
for fixed 4>h and E\, the electric field strength at the 
emitter surface Eq has a square root dependence on the 
small quantity h c — h. For h = h c , the electric field at 
the emitter surface is completely screened by the space 
charge. At given </>& and E\, the space charge density 
decreases as the interelectrode distance h is reduced. The 
opposite limit, where the space charge is absent (Eq = E\ 
and 4>{z) = E\z), is reached at h = (3/4)/i c . 

So, in the framework of the discussed analogy we can 
identify (i) the angles 9 m , 9 C , and 9t with the interelec- 
trode distances h, h c , and hx = (3/4)/i C) respectively, 
(ii) the electric field strength on the cone surface at a 
unit distance from the apex E with the strength at the 
surface of the planar emitter surface Eq, (iii) the electric 
field strength at the symmetry axis at a unit distance 
from the cone apex |V</>| = 1/2 with the strength at the 
opposite electrode E\, and also (iv) unit potential dif- 
ference between the cone axis and its surface for r = 1 
with the fixed potential difference <ph for a plane vacuum 
diode. In order to compare the characteristics of conical 
and planar diodes, we apply the linear mapping of the 
angle interval 9x < 9 m < 9 C into the interval of distances 
hx < h < h c , 

h~h T _ 9 m - 9 C 
h — h c 9 m — 9t 

The analogs of the functions F(a) and E(a) correspond- 
ing to the one-dimensional model (|2"5)l , 

F pI = F ( a ^ > E pI = IT E M > ( 33 ) 

are presented in Fig. [3J It can be seen that the properties 
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FIG. 5: The current J as a function of the applied potential 
difference V. The dotted line shows the asymptotic shape of 
the current- voltage characteristic corresponding to the Child- 
Langmuir law. 

of conical (fF2|) and planar (|28|) models are qualitatively 
similar. Consequently, we should expect that E ~ (9 C — 
9m) 1 / 2 at 9 m — > 9 C , or, in terms of the cone half-angle, 
E ~ (a — Uc) 1 ! 2 at a — > a c , as we wished to show. 

Taking into account such a dependence near the lim- 
iting regime (|25|) . we approximate the calculated depen- 
dence of E on a, by the formula 

E(a) = -0.9754£ - 0.6176C 2 + 0.8342C 3 - 0.2164£ 4 , 

(34) 

1 /2 

where £ = [(a — a c ) / (ar — ot c )] ■ The error of the 
approximation (|3~4"]) is less than 0.1% over the range 
a c < a < olt- According to the balance condition (fTl)]) . 
the obtained dependence E(a) determines the relation 
between geometry of the emitting cusp, i.e., the angle a, 
and the external control parameter of the system, i.e., 
the potential V ~ <J>l . 

D. The relations between model parameters 

The quantity F(a), which is determined by (Tl9l) . spec- 
ifies a particle flux from the top part of the cone, < r < 
1. The angular dependence of this function obtained as 
a result of numerical solution of the problem (fT2)) -(IT4 )) is 
presented in Fig. El The calculated function F(a) can be 
approximated by the relation 

F(a) ~ 1.8863C + 0.0956C 2 - 0.2259C 3 - 0.1167C 4 , (35) 

where £(a) = 1 — £ 2 = (ax — a) /(ax ~ a c ). The error 
of this approximation is less than 0.15%. Note that the 
function F(a) monotonically decreases from ~ 1.637 to 
zero as the angle a changes from a c to cut- 
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FIG. 6: The change in the cone half-angle versus the emission current (left) and the current-voltage characteristics (right) for 
a gallium liquid-metal ion source. The points correspond to the experimental data of and the lines to the theory. 



The expressions (fT6|) and ([T9|) together with (|34|) and 
(f3"5"j) allow us to obtain the dependence of the cone half- 
angle a on the current J as well as the current-voltage 
characteristic of the cone, i.e., the dependence of J on V. 
These relationships are plotted in Figs. |4] and [5j respec- 
tively. 

From Fig. |4] it can be seen that the angle monoton- 
ically decreases with increasing current. It is equal to 
the Taylor angle or ~ 49.29° at zero current and tends 
to the angle a c ps 21.89° in the formal limit of an infi- 
nite current. Note that the interpretation of experiments 
in the framework of our model corresponds to 
the dimensionless current in the range < J < 1 (see 
the inset of Fig. [4]). For larger current values, the cone 
structure becomes instable. 

An important feature of the theoretical current-voltage 
characteristic (Fig. [5]) is its threshold character. There is 
no electric current, J = 0, if V < Vo — y/ctgar / f (ar) — 
0.9512. This is related to the impossibility of a balance 
between the electrostatic and capillary forces at a rel- 
atively small potential difference. The capillary force 
will dominate and the cone structure will break. In the 
formal limit of large V the model yields the universal 
Child-Langmuir law: J — > F(a c )V 3 ^ 2 , which describes 
the regime of space-charge limitation of the current as 
a result of complete screening of the electric field at the 
electrode surface. 

Excluding the parameter a from (|T6| and (IT9l , we can 
rewrite the dependence of J on V in the explicit form. 
The function F linearly goes up to the value F(a c ) as 
a — » a c (see Fig. [3]), and, as discussed above, the function 
E tends to zero according to the square root law. As a 
consequence, the dependence of J on V must have the 
form J w joV 3 / 2 + jiV^ 1 / 2 in the limit V -} oo. With 
this estimate taken into account, for V > Vq the required 



dependence can be approximated by the expression 



J2_ 

1/4 



21 



(36) 



The coefficients of the approximation (|36p are the fol- 
lowing: j = 1.6372, ji = -1.8635, j 2 = 0.3523, 
73 = —0.00676; the approximation error is less than 
0.00015 in the absolute value. 

Thus, the relations (|16l) and p^|) together with the ap- 
proximate expressions (l34l) . (|35|) . and (l36l) completely de- 
termine the integral characteristics of the system. They 
establish relationships between the following quantities: 
the applied potential difference V, the emission current 
J, and the cone half-angle a. This will allow us to com- 
pare our theoretical calculations with available experi- 
mental data in the next Section. 



III. COMPARISON WITH EXPERIMENTAL 
DATA 

In order to compare the results of our calculations with 
experimental data, we assume that the control parameter 
of the model $l is directly proportional to the potential 
difference U applied to the interelectrode space, i.e., 



kU, 



(37) 



where the proportionality coefficient k is determined by 
constructional features of an experimental facility and 
does not depend on the emission current. Then, after 
switching to dimensional quantities, the expressions (| 16[) 
and (|19[) become 



^. '2a\^^ 2 yjL(a) cot a 
so J nE{a) 



(38) 
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FIG. 7: The change in the cone half-angle versus the emission current (left) and the current-voltage characteristics (right) for 
indium and tin liquid-metal sources. The points correspond to the experimental data of [15 . [13] and the lines to the theory; 
the solid line corresponds to indium, the dashed line to tin, and the dotted line to tin under the condition of the emission of 
singly-charged ions Sir + ' only. 



I = £ 



2q 



1/2 



F(a) {kU) 3/2 . 



(39) 



Here we take into account that the characteristic size L 
of the top part of the infinite model cone (the electric 
current from this part of the cone is identified with the 
emission current from the experimentally observed coni- 
cal spike) can depend on the value of the emission current 
and, as a consequence, on the angle a. 

The relations (f38]) and (|39|) together with the approx- 
imate formulas (|3"5|) or (|3"6")l allow one to de- 
termine the free parameters of the model, namely, the 
characteristic size of the cone L(a) and the coefficient k 
characterizing the electric field distribution in the exper- 
imental facilities, from the available experimental data 
on current-voltage and current-angle dependencies for 
liquid-metal ion sources. We have used data for emis- 
sion into vacuum from liquid indium [Til ], gall ium [l5j . 
tin and also from liquid alloys Au+Si [16j], Co+Gc 
[H|, and Au+Ge [T^. In our calculations we have taken 
the following values of the surface tension: a = 0.572 
Im (In, pj), 0.735 N/m (Ga, 0), 0.560 N/m (Sn, 

1.84 N/m (Co+Ge, pi]), 
The angular dependence 



), 2.20 N/m (Au+Si, 
and 1.62 N/m (Au+Gc, 
of L was approximated by the two-parameter function 



19] 
19) 



L(a) = L C + (L T - L c ) 



CUT 



(40) 



where Lt and L c are parameters. It should be noted 
that the maximum half-angle a max , corresponding to zero 
emission current in the experiments [14| - tl9| . slightly dif- 
fers from the Taylor angle cut- Then, for comparison with 
the theory, the experimentally observed values of a have 
been corrected by ar — cc max . 



Due to its physical properties (low temperature of 
melting, low pressure of saturated vapor, primarily single 
ionization in the emission processes), gallium is the most 
convenient and reliable metal for investigating the oper- 
ation of liquid-metal ion sources Q . In the work [l5| the 
gallium ion sources, including a tungsten needles with tip 
radius Rt — 2 /im and R t — 20 /im, were investigated. As 
a result of treating data on the current-voltage charac- 
teristics and geometry of the observed structures on the 
surface of liquid gallium, presented in [l5j . we obtain: 
k = 0.0905, L T = 0.65 /mi, L c = 0.84 /tm (R t = 2 /mi); 
k = 0.076, Lt = 0.84 fxm, L c = 0.84 /im (R t = 20 /mi). 
The achieved quality of describing the experimental data 
is demonstrated in Fig. [B] Application of a sharper nee- 
dle (with R t — 2 /im) provides a stronger focusing of 
the electric field that corresponds to a larger value of the 
coefficient k and, consequently, to a smaller value of the 
threshold potential difference (see Fig. EJ)- On the other 
hand, a small initial radius of curvature probably begins 
to "constrict" the developing cone; the value of Lt cor- 
responding to R t — 2 /jm is somewhat smaller than for 
Rt = 20 /an. According to the relations (|3"5|) and (|3"9"1) . 
smaller size of the cone leads to faster change of the cone 
half- angle with increase in emission current, all other pa- 
rameters (the surface tension coefficient, the mass and 
charge of emitted particles) being the same (see Fig. [S]). 
Note that the obtained values of Lt are comparable with 
the size of the experimentally observed cones (they are 
of the order of 1 /tm) . 

Fig. [7] shows the current-voltage and current-angle de- 
pendencies for the conical spikes observed at the surface 
of liquid indium [Tij and tin 17 1. Probably, the strongly 



nonlinear character of the current-voltage dependence for 
indium is an experimental error [Tij . In later publications 
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FIG. 8: The change in the cone half-angle versus the emission current (left) and the current-voltage characteristics (right) for 
Au+Si [lfj], Co+Ge [l8|. and Au+Ge [l^] alloy liquid-metal ion sources. The points correspond to the experimental data and 
the lines to the theory. 



[18l Il9j the same authors approximate this dependence 
by a linear function; its slope is approximately the same 
as for the corresponding theoretical curve in Fig. [7J 

It has been noted in [17j that both singly-charged ions 
Sn( + ) and doubly-charged ions Sn( 2+ ) are emitted from 
the surface of tin. The fraction of doubly-charged ions 
is x 2 = N 2 /{N l + N 2 ) ~ 45% HU. It is easy to verify 
that the relations obtained above are applicable for a 
flow consisting of k different types of particles. The only 
necessary modification refers to the expression for total 
current (j3"9")l . We should apply the following change: 




E 



(41) 



where qi, rrii, and Xi are the charge, mass, and relative 
fraction of particles of i-type. So, for the case of tin 
(<72 = 2gi = 2e, m 2 = m\ — m), the expression (|39[| 
transforms into 



1 + (2y/2 -l)x 2 /2e\ 1/2 3/2 

I = — £ — F(a) (kU) . 

1 + x 2 \m I 



(42) 



where e is the elementary charge. Indium ions are emit- 
ted in the singly- charged state In( + ) [l9j], so that such 
modification is not required for an indium liquid-metal 
ion source. The theoretical curves presented in Fig. [7] 
correspond to the following values of the parameters: 
k = 0.0495, L T = 0.4 /itm, L c = 0.3 fan (In); « = 0.075, 
Lt — 1.1 (im, L c = 0.78 /im (Sn). 

Indium and tin have close atomic weights and surface 
tension coefficients. In this connection, a sharp distinc- 
tion between the dependences I(U) and a(I) correspond- 
ing to these metals is surprising 17]. One of the likely 
reasons of the different behavior of indium and tin is the 



presence of a large amount of ions Sn 2+ in the experi- 
ments [171 ]. Ions with larger charge numbers move faster 
away from the electrode and, as a consequence, provide 
less screening effect. Therefore, an increase in the voltage 
is accompanied by a larger increase in the emission cur- 
rent and by a relatively slow change in the cone shape. 
The dependences I(U) and a(I), corresponding to the 
assumptions that only the singly-charged ions Sn( + ) are 
emitted from liquid tin and the parameters k, Lt, and L c 
are fixed, are shown in Fig.[7]by the dotted line. One can 
see that considering only singly-charged ions slightly ap- 
proaches the model characteristics of liquid-metal emit- 
ters based on liquid tin and indium. However, the initial 
difference between the curves is too large and cannot be 
explained only by the presence of ions Sn^ 2+ ) . 

Another reason of the behavior difference between 
liquid-metal ion sources is in the conditions of experimen- 
tal studies [lij and [l?} • The theoretical model presented 
in this work points to the essential difference between the 
characteristic sizes Lt of the cones. One can assume that 
the developing cone was "constricted" by the small size 
of the tip of a tungsten needle in the experiment with 
indium (i? t = 1.0 fan according to [l4| )• Although the 
value of the tip radius Rt = 4 fim for the experiments 
with indium was indicated in the subsequent paper [l9| . 
it was measured with the help of the image of the liq- 
uid surface instead of the needle. The above analysis 
of gallium emitters with different needles shows that the 
"constriction" effect becomes appreciable if R t < 2 fan. 

Finally, Fig. [5] shows the dependencies I(U) and a(I) 
for liquid alloys Au+Si [16], Co+Ge [T|, and Au+Ge 
that correspond to the model coefficients k — 0.084, 
L T = 1.0 fan, L c = 0.45 /im (Au+Si); k = 0.084, 
L T = 0.27 fan, L c = 0.08 ^m (Co+Ge); n = 0.051, 
L T = 0.30 fim, L c = 0.10 fim (Au+Ge). Theory and ob- 
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servations are in rather good agreement only for Au+Si 
from the above-listed alloys. The alloys containing ger- 
manium are distinguished by the presence of breaks on 
the experimental curves a(I). Such a behavior cannot 
be described in the framework of the proposed model. 
The anomalous character of the dependence a (J) can be 
caused by changing properties of charged particles flow. 
So, "freezing" of the cone angle along with the continu- 
ing growth of the emission current, which is characteristic 
for the mentioned alloys (see Fig. [5J , allows us to sup- 
pose that, from some point, the average ratio of mass 
and charge dencities over the flow essentially decreases. 
Unfortunately, the experimental works [HI, [l!| does not 
include information about changes in the mass/charge 
ratio with the emission current. In Fig. |8] the theoretical 
curves are obtained under the assumption that the inten- 
sive emission of doubly-charged germanium ions Ge( 2+ ) 
starts with some value of the applied potential difference, 
U\ (or current, The quantities U%, I± are of the or- 
der of 3.87 kV, 30 mkA for the alloy Co+Ge, and 6.3 kV, 
15 mkA for Au+Ge. In both cases, the potential drop is 
<I>l = kU ~ 0.32 kV at the distance of the order of the 
cone size (R = L). 



been obtained. This allows us to find the dependences 
of the cone angle and emission current on the applied 
potential difference. They correctly represent the main 
features of the operation of a liquid-metal ion source. A 
comparison was made between the developed theoretical 
model and available experimental data for emission from 
pure gallium, indium, tin, or from alloys Au+Si, Co+Ge, 
and Au+Ge. Basing on our theoretical results, we have 
proposed explanations for some specific features of the 
emissive behavior of these systems. So, the difference in 
behavior between indium and tin having close character- 
istics (the coefficient of surface tension, the mass/charge 
ratio) is probably related with the presence of doubly- 
charged ions Sn' 2+ ' for tin emission and with use of a 
sharper tungsten needle in the experiments with indium. 

In conclusion, let us note once again that the devel- 
oped model does not pretend to describe all aspects of 
operation of liquid-metal ion sources. In particular, the 
consideration of the jet-like protrusion at the cone ver- 
tex, from where the vast majority of ions are evaporated, 
remains beyond the scope of the model. The self-similar 
solutions used above describe LMIS operation in terms 
of the averaged characteristics (a, /, and U) only. 



Conclusion 

In the present paper we have developed the self- 
consistent model describing how the space charge near 
the emitting cone apex affects its shape. Our approach 
is based on the self-similar reduction of the equations, 
which govern the spatial distributions of the electric field, 
ion velocity field, and particle concentration, to the sys- 
tem of ordinary differential equations. As a result of 
numerical solution of this system, the conditions of the 
mutual compensation of the capillary and electrostatic 
forces on the conic surface of a liquid-metal anode have 



Acknowledgments 

The study was performed within the framework of the 
Program of Interdisciplinary Projects between the Ural 
Division and Siberian Branch of the Russian Academy of 
Sciences and the Program "Basic Problems of Nonlinear 
Dynamics" of the Presidium of the Russian Academy of 
Sciences. It was financially supported by the Russian 
Foundation for Basic Research (project 07-02-96035), by 
the Presidential Program of Grants in Science (project 
MD-2553. 2007.2), and by the Dynasty Foundation. 



[1] J. Zeleny, Phys. Rev. 10, 1 (1917). 
[2] G. I. Taylor, Proc. R. Soc. A 280, 383 (1964). 
[3] R. G. Forbes, Vacuum 48, 85 (1997). 
[4] L. Oddershede and S. R. Nagel, Phys. Rev. Lett. 85, 1234 
(2000). 

[5] N. M. Zubarev, JETP Lett. 73, 544 (2001). 

[6] N. M. Zubarev, Phys. Rev. E 65, 055301(R) (2002). 

[7] V. G. Suvorov and N. M. Zubarev, J. Phys. D: Appl. 

Phys. 37, 289 (2004). 
[8] R. Cornier, Appl. Phys. 19, 365 (1979). 
[9] N. K. Kang and L. W. Swanson, Appl. Phys. A 30, 95 

(1983). 

[10] D. R. Kingham and L. W. Swanson, Appl. Phys. A 34, 
123 (1984). 

[11] G. S. Cho, Y. Seo, and S. O. Kang, J. Appl. Phys. 72, 
5892 (1992). 

[12] P. Wang, S. Maheshwari, and H.-C. Chang, Phys. Rev. 

Lett. 96, 254502 (2006). 
[13] G. I. Taylor, Proc. R. Soc. London 291, 145 (1966). 



[14] B. Praprotnik, W. Driesel, C. Deitzsch, and H. Neidrig, 

Surf. Sci. 314, 353 (1994). 
[15] W. Driesel, C. Dietzsch, H. Niedrig, and B. Praprotnik, 

Ultramicroscopy 57, 45 (1995). 
[16] W. Driesel, C. Dietzsch, and R. Miihle, J. Phys. D: Appl. 

Phys. 28, 787 (1995). 
[17] W. Driesel and C. Dietzsch, Appl. Surf. Sci. 93, 179 

(1996). 

[18] W. Driesel, C. Dietzsch, E. Hesse, L. Bischoff, and J. Te- 
ichert, J. Vac. Sci. Technol. B 14, 1621 (1996). 

[19] W. Driesel, C. Dietzsch, and R. Miihle, J. Vac. Sci. Tech- 
nol. B 14, 3367 (1996). 

[20] C. D. Child, Phys. Rev. 32, 492 (1911). 

[21] I. Langmuir, Phys. Rev. 2, 450 (1913). 

[22] G. L. R. Mair, J. Phys. D: Appl. Phys. 17, 2323 (1984). 

[23] G. L. R. Mair and R. G. Forbes, J. Phys. D: Appl. Phys. 
24, 2217 (1991). 

[24] G. L. R. Mair and R. G. Forbes, Surf. Sci. 266, 180 
(1992). 



11 



[25] P. T. Kirstein, G. S. Kino, and W. E. Waters, Space- 
Charge Flow (McGraw-Hill, New York, 1967). 

[26] J. M. Finn, T. M. Antonsen, and W. M. Manheimer, 
IEEE Trans. Plasma Sci. 16, 281 (1988). 

[27] G. S. Boltachev and N. M. Zubarev, Europhys. Lett. 76, 
36 (2006). 

[28] L. Landau and E. Lifshitz, Course of Theoretical Physics, 
vol. 6: Fluid Mechanics (Pergamon, Oxford, 1987). 



[29] A. M. Ganan-Calvo, Phys. Rev. L 79, 217 (1997). 

[30] A. Barrero, A. M. Ganan-Calvo, J. Davila, A. Palacio, 
and E. Gomez-Gonzalez, Phys. Rev. E 58, 7309 (1998). 

[31] A. Hirabayashi and J. F. de la Mora, International Jour- 
nal of Mass Spectrometry and Ion Processes 175, 277 
(1998). 

[32] E. W. Muller, Phys. Rev. 102, 618 (1956). 



